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Abstract 

A graded Holstein-Primakoff realization of the SU(2/1) algebra is pro- 
posed. A spin-wave theory with a condition that the sublattice magnetization 
is zero is discussed. It is shown that the generalized spin-wave theory is ap- 
propriate to investigate the long-range spiral (incommensurate) phase of the 
t — J model at r = 0. The spin-spin correlator is calculated. 
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I. Introduction 

Among the many electronic models which are being currently studied in 
the context of high-Tc superconductivity, the two dimensional t — J model 
is the simplest one and it is rather general to describe the basic physics of 
the new superconductor^. It is an effective model for the large U Hubbard 
model with J = and there are arguments that one can derive it from a 
multiband Hubbard model^'^. 

At half-filling the t — J model reduces to the Heisenberg model. This is a 
welcome result, since, the methods developed for the Heisenberg model can be 
extended and applied to investigate the properties of the t — J model, atleast 
at nearly half-filling. Schwinger boson slave-fermion mean field theory was 
used and an antiferromagnetic ordered Neel state at half-filling was found^. 
This state evolves into long-range spiral (incommensurate) antiferromagnetic 
states at large doping. Quantitatively, this is described by the dependence of 
the spirality angle Q on the doping^. Originally, the spiral phase of the t — J 
model was obtained by Shraiman and Siggia who used a somewhat different 
approach^. 

Recently, Takahashi has formulated a modified spin- wave theory of Heisen- 
berg (anti)ferromagnets'''^. He has supplemented the usual spin- wave theory 
with the constraint that the magnetization of each site is zero. This enforces 
the condition that the total number, on average, of spin waves per site is 5, 
and that the sublattice rotational symmetry is not broken. In one dimension 
the modified spin wave theory yields excellent agreement with Bethe ansatz 
results for spin S = 1/2. Takahashi's results are in quantitative agreement 
with the Schwinger boson mean field theory^ and with the renormalization 
group theory^. 

In this paper we formulate, along the same line, a generalized spin wave 
theory appropriate to investigate the spiral phase oft— J model. Some results 
of Ref. 10 are used. 

In Section II the canonical bose- and fermi- operators are used to realize 
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the SU(2/1) algebra. The graded Holstein-Primakoff representation is writ- 
ten in a local frame. The effective hamiltonian is obtained in leading order 
of S~^, as a function of the angles 6i which determine the local frame. To 
preserve the sublattice rotational symmetry we impose two additional con- 
ditions. To enforce them, we introduce two terms in the hamiltonian with 
chemical potentials A and For the special values oi di = kir, ( is equal 
to zero and we reobtain the Takahashi's spin-wave theory. We express the 
angles 6i by the spirality angle Q. To determine the last one, we calculate 
the energy of the system per site. The physical value of Q is that which 
minimizes the energy. It depends on the doping. 

In Section III we calculate the large distance asymptotic of the spin- 
spin correlation function. At T = we obtain long-range correlation which 
is characterized by the angle Q and the quantity m which is related to the 
spin-wave bose-condensate. At finite temperature, the correlation function 
exponentially falls when the separation between the two sites is large. At 
small doping the spin-spin correlator looks like in the theory of Heisenberg 
antiferromagnet. When the doping is large enough it looks like in the theory 
of Heisenberg ferromagnet. 

Section IV is devoted to the concluding remarks. 

II. Graded spin-wave theory 

The t — J model is defined by the Hamiltonian 



where c^(cjo-) are the fermi operators of the electron {a = 1, 2) on site 

of a 2D square lattice which act on states with no double occupancy on a 

— # 

lattice site Si are spin operators, and fi is the chemical potential. By 
<i,j> we denote the sum over the nearest neighbors. 

The eight operators Qo-, c^, S^, Sf and form a basis of the SU(2/1) 
graded algebra. It is most convenient to rewrite them in terms of the Hubbard 



h = -t y: [ 



<i,j> 



cl^Cja + h.c] + J [Si - Sj- -UiTlj 

<i,j> \ 4 




(1) 
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operators X^^ {a,b = 1,2, 3) 



= Xf^ 



a = 1,2 



Tin 



— A- +A- ^- — 2 (^A ■ -A, 
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C+ _ vl2 C- _ v21 

Oj — — 

They satisfy the following graded commutation rules 



(2) 



(3) 



The graded algebra SU(2/1) can be realized using a canonical bose a^, 
and fermi tpi, ip^' operators 



CLj 



= 6ij; [ipi,ipf]_^ = 6ij 



(4) 



Cil = 



(2S - a^ai - ip^ipi) 



1/2 0. 



COS — h flj sm — 
2 2 



Ci2 = 



^2 = 



(25 - a.+ttj - V'i^V'i)^^^ cos ^ + sin ^ 

1/2 . 9. 



1/2 . 9 



1 9i 
sin \- ai cos — 

2 2 



(2S — ajai — ij), 



sin — + at cos — 
2 * 2 



= 25 - 



1 



St = - sin 9i 



(25 — a^flj — i/ji'ipiY^'^ ai + (25 — a^ai — ip^'ipiy^'^ 
+ cos 9i (^S - aftti - ^i/^fipi^ 
Si' = cos^ ^ (25 — a^ai — i/j/ipiY^'^ — sin^ (25 — a^Oj — V'l^V'i) 

- sin 6*^ 1^5 - a^ai - ^V'i^V'i) 



(5) 



1/2 
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Si = cos^ (2S — afai — iplipiY^^ — sin^ ^ (2S — afai — 



- sin 6'i ( 5 - alai — -il^^i/Ji 
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where 6i are arbitrary angles which run over intervals with length tt. When 
6i = 0,7r, 27r... the Eqs. (5) reduce to the graded Holstein-Primakoff repre- 
sentation^^. It is easy to prove that so defined operators satisfy the algebra 
Eqs. (2,3). In particular, the spin operators satisfy the SU(2) algebra and S 
parametrizes the representation of the SU(2/1) algebra. At half-filling S is 
exactly the spin of the system. The physical relevant case is 5* = |. 

1 /2 

In spin-wave approximation one must replace (2S — a^ai — '(l^t'^i) by 
{2Sfl'^ [l-js {alai + ipti^i)] in Eqs. (5). Then, in leading order of S''^, the 
Hamiltonian Eq. (1) takes the form 



where 



and 



h = hce-\- hq (6) 



hce = JS'' E C0S(^, - Oj) (7) 
<i,j> 



hq = —JS X] cos(6'j — 6j){aiai + a-aj)^ 

<i,j> 



JS 

E [cos(6'i - 9j) + 1] {a^aj + ajai) 



+ — E [cos{ei-ej)-l]{aiaj^aial) (8) 



2 

JS 

2 <i,j> 
JQ 

E [cos(^. - e,) - 1] {^lJt^^^^ + V^;^^,) 

^ <i,3> 

+2St E COS ^^(V'+V'i + ^t^i) -f^U^S - ^t^i) 

<i,j> ^ i 

Farther on, we shall use special values for the angles 9i 

Oi = Q-n; Q^ = Qy = Q (9) 
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where 'k <Q <2ti and = (rf , rf) are the coordinates of the lattice site i. 

The Eqs. (6-8) yield the conventional spin wave theory. To preserve 
the sublattice rotational symmetry one must impose, by hand, an additional 
condition. Here it is 

<ni-Si>=Q (10) 
where Hi is a unit vector which fixes the local frame we have chosen in Eqs. (5). 

ni = (— sin^j, 0, cos^j) (11) 

When 9i = 2/c7r the Eq.(lO) is the condition imposed by Takahashi for 
ferromagnets^. When Oi = kn Eq. (10) coincides with the condition dis- 
cussed in spin-wave theory of antiferromagnets^'^^ 

Making use of Eq. (11) and the representation for the spin operators 
Eqs. (5), one can rewrite the condition Eq. (10) in the form 

< a+ai > +^ < '^t'ipi >= S (12) 

To enforce the constraint Eq. (12) we introduce a new term in the Hamilto- 
nian 

h^h-XYl {2a+ai + i/j+iPi - 2S) (13) 

i 

with Lagrange multiplier to be determined by Eq. (12). 

It is not difficult to check, that Eq. (12) has a solution just for Q = ir, 27r. 
This is easy to be understood because the identities 

< ttitti >=< a^al >= (14) 

are satisfied only for these values oi Q. To ensure the implementation of the 
Eq. (12) for arbitrary Q in the interval [tt, 27r], one must impose in addition 
the condition Eq. (14) also. We introduce in the Hamiltonian a term with 
Lagrange multiplier A which enforces the condition Eq.(12) and a second term 
with Lagrange multiplier C, which enforces the conditions Eq. (14). The full 
Hamiltonian in the momentum space reads 

h = 2JS'^N cos Q + 2SN{X - //) (15) 
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where 



k 



Cbik) = -2A - 45 J cos Q + 25 J(l + cos Q)jk 
Dk = SJ[cosQ-l]jk + C 
ef{k) = II- X + 25J(1 - cos Q) + 8St cos ^^7^ 



7fc = -(cos /Ca; + COS (16) 



1 
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N is the number of lattice sites, and k runs over the Brillouin zone. Through- 
out this paper we put the lattice spacing equal to one. 

The Bose part of the Hamiltonian Eq. (15) is diagonalized by a Bogoli 
ubov transformation. The result is 

/i6 = E \Eb{k)alak + Eo{k)\ 

k 

where 

^21 1/2 



E^ik) = [el{k) - ADl 

Eo{k) = \ [E,{k) - e,{k)] (17) 
The free energy of the system is given by the expression 

F = 2 J5^ cos Q + 25(A - /i) + ^ E Eo{k) (18) 

^ k 

i- E In (1 - e-^^^^'^)) - -1- ^ In (1 + e'^^/^^)) 



where (5 = (ksT)'^ is the inverse temperature. The three chemical potentials 
A, C, and /i are determined by the equations 

dF dF dF 

ax='' 

It is more convenient to write their linear combinations 



TV 5 1 + e-/^^/W ^ ^ ^^^^ 
where (5 is the doping parameter. 

We introduce new chemical potentials 771, 772 which are determined by 
the equalities 

2SJr]i = -A-C-25JcosQ 

25 J772 = - A + C - 25 J cos Q (22) 



Then 



e6(A;) - 2Dk = 4SJ{r]i + 7^) 
Cbik) + 2i:»fc = 45 J(772 + cos Q-fk) (23) 



£;j,(A;) = 45JV(r/i + 7fc)(^2 + cosQjk) 

From Eqs. (19-23) it follows, that when Q = ir, 27r, the solution is rji = r]2 
or ( = 0. This means, that for these special values of Q the identifies Eq. 
(14) are satisfied without introducing a second term in the Hamiltonian and 
one reobtains the Takahashi's spin-wave theory. 

The chemical potentials depend on the temperature and rjiiT) > 1; 
^2(7") >| cosQ |. Next, we analyze the low-temperature properties of Eqs. 
(19,20) r]i{0) and r]2{0) depend on Q and 6. For given doping 6 we shall 
consider three cases. 

When IT < Q < Qi < ^TT, we obtain from Eqs. (19,20) (see Takahashi^) 

25+l-(5-1.394 I cosg 1^/^ + l ^^ ln2{4SJ(3f \ cosQ | (771-I) = (24) 

zttdJ/j 

25+l-(5-1.394 I cosQ T^/^ +--— — In 2(45 J/?) ' (772- | cosQ |) = 

ZnbJp I cosQ I 

From Eq. (24) one gets 

_ 1 = I -47r5Jmi/3 

" 2i4Sjpy\cosQ\ 



772- I COS (5 1 = 



1 



2(45J/3)2 



3 — Aw S J 1712 1 COS Q I /? 



1 1/2 
-1/2 



where 

2mi = 25 + 1 - (5 - 1.394 | cosQ 

2m2 = 25 + 1 - J - 1.394 | cosQ | 

mi increases and 1712 decreases as | cosQ \ decreases, and 1712 
upper bound of the interval 

/ 1.394 ^21 



Q{i) 



arccos 



.25 + 1 -(5. 



When ^ < (^2 < Q < 27r one obtains form Eqs. (19,20) 



(25) 



(26) 
for the 



2S+1-6-{cosQY/^ = 



1 



27rSJ(3 



In 185 JAi 



(Ai + A2)V^ + (Ai + 4A2)V^ 2 

[2(Ai + A2)V2 + (Ai + 4A2)V2J (cosg)V2 



1 



27rS J cos Q/3 



where 



25+l-J-(cosQ)-V2_ 
ln25J[(Ai + A2)i/2 + (Ai + 4A2)^/2] 

Ai = {rji - l){r]2 - cos (5) cos (5 
A2 = ^{m - cosQrjif 



(cosQ)V2 



(27) 



(28) 



Straightforward calculations lead to the following low-temperature assymp- 
totics of the chemical potentials r]i and 772 

-(Mi-M2)/3\ 1/2" 



771 - 1 = 



772 - COS (5 = 



25J(cosg)V2/3 

(cosQ)i/2e-^^^ 



25J/3 



3 ^ g-(Mi-M2)/3 ^i_g- 

9 _ g-(Mi-M2)/3 

3 + g-(Mi-M2)/3 



9 _ g-(Mi-M2)/3 



+ 



1 



-(Mi-M2)/3\ 1/2- 



9 _ g-(Mi-M2)/3 ' I 9 _ g-(Mi-M2)/3 



(29) 
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where 

Ml = AttSJuii; M2 = ^'nSJm2 cos Q (30) 
and mi, 777-2 are given by 

2mi = 2S+l-5-{cosQf''^ 

2m2 = 25 + 1 - (5 - (cos g)-^/^ (31) 
For Q = 27r, ttii = ?t72 = ra^ Mi = M2, rji = r]2 = rj and 

For the other values of Q, M2 < Mi. Q2 is defined from Eq. (31) with 777,2 = 

Finally, let Q runs the interval Qi < Q < Q2- Then 771 (T) goes to one, 
when T goes to zero, and 772(0) >| cosQ \. 

[AS J mm- cos qf 



where 7711 and 772(0) are determined by the equations 

1/2 

7fc 



2S + 1-S = ^Y.( ] (35) 

N ^ \t]2 ^ cos Q-fk J 

and 7772 = 0. 

Let us recapitulate the results we have received above. When the tem- 
perature goes to zero 771 (T) goes to one for all values of Q. At the same time 
7/2(0) =1 cosQ I if TT < Q < Qi, and 7771,7772 can be found from Eqs.(26). 
When Qi < Q < (^2, ^2 = and 7/2(0) and 7?7i can be obtained from the 
system Eqs.(35). Finally, when Q2 < Q < 27r, 772 = cosQ and 7771,7772 are 
given by Eqs. (31). 
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For the further convenience we have used a common notation for the 
variables mi, 1712 when ir < Q < Qi and when Q2 < Q < 27r, but their 
origin is different. At zero temperature, when 'n < Q < Qi, r]i = \ and 
r}2 =1 cosQ I, the bose system condenses at the wave vectors k* = (±7r, ±7r) 
and kg = (0,0). Then, the spin- wave correlators read 

1 



< ak >= - 



1 



(1 - h,k*) (1 - h.ko) 



+ ■^'nk*Ok,k* + -jnk„Ok,ko 



(36) 



< a^a-k >=< CLk^-k 
•-^|ycth^E6(A;)(l - 6k,k*){X - h,ko) - ]^nk*5k,k* + ^ntjkx 



and 



1 ^ 

= irrl^^k*: 
iV k* 



1712 



N 



(37) 



When Q2 < Q < 27r, ryi = 1 and r]2 = cosQ the bose system condenses at 
the wave-vector k* and for the spin-wave correlators are obtains 



+ 1 



eb{k) 3 



1 



.1 — h,k*) + -;^nk*^k,k* 



< aktt-k >=< CLkaZk > 
where 



•-^|ycth^E6(A;) (1 - dk.k*) - ^nf*Sk,k* 



(38) 



1 ^/ (1) , (2)x 



2N 



-^l^K^k* -^k*) - ^2 



(39) 



When Qi < Q < Q2 the spin- wave correlators are given by Eq.(38) with 
n<l' = ni? {m, = 0). 

To complete our analysis we shall discuss the magnon spectrum as a 
function of Q. From Eqs (23) we obtain that when ir < Q < Qi there are 
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two antiferromagnetic magnons: Ei,{k) c::^ Vs \ k — k* \ when | k — k* |^ 
and Ef){k) Vs \ k — k'^ \ when | A; — A;^ |— ^ 0, with spin-wave velocity 
Vs = 4:SJy/\ cosQ I /2. When Q exceeds Qi one of the magnons opens a 
gap Eb ~ Vs\J \ k — k^ \^ where = 4(r/2— | cosQ |)/ | cosQ |. In 
the interval (37r/2,Q2) there is only one magnon E}j c^l v'^ \ k — k* \ with 
v'g = 2SJ{r]2 — cosQ)^/^. The spin- wave velocity v'g approaches zero when 
Q approaches Q2 and in the interval [Q2, 27r] we obtain the ferromagnetic 
magnon E^ ~ SJ^/cosQ \ k — k* p. 

The angle Q is a free parameter in the theory. To determine it we shall 
calculate the energy of the system per site. 

From Eqs.(5), taking into account the conditions Eqs. (12,14), one ob- 
tains in spin-wave approximation 



cos{di—dj) 



< Si ■ Sj >= 



< ajaj >< ttj'tti > + < a^ttj' >< a^a^ > 



1 

4 



(40) 



and 



— cos 



9i — 6j 



< C+ Cj^ + /i.e. > = 

{S + 1-6) (< ij+il^i > + < ^P+i/jj >) + 



< niUj >= {2S - 5f- < ijjjiljj >< i/j^i/ji > 



(41) 



Putting Eqs. (40-42) into Eq.(l) one obtains for the energy per site at 
T = and TT < Q < Q2 



E = -^(25 -5)2 + 1(1- cos Q) 



^E7.n-e/(A:)) 



+ 



J cos Q 



2s+i-^-li:(i-,,)( 



1/2- 



772 + COS Qjk. 



+ 



12 



J cosQ 



2S 



r]2 + cos QjkV^^ 



+ 7fc 



Q 

+2t cos — 
2 



25^i:(i-7.)( 



+ 

1/2 



2iV 
l + 7fc 



1/2 



7^2 + COsQjk. 



+ 25 + (5 



(42) 



where 772 is defined above and 0{—ef) is the step function. 
When Q2 ^ Q ^ 27r, the energy is given by 



E = -^(25 - 6)^ + hi - cosQ) 



^E7fc^(-e/(fc)) 



J cos Q 
4 

JcosQ 



25 + 1 
25 + 1 



+2i cos 



^E7fc^(-e/(/c)) 



5 



26 



1 



Vcos(5 



cos Q 
1 + 



+ 



+ 



1 



2a/cos Q 



+ 



VcosQ 



For a given 6 E is a, function of Q, and the physical value of Q{Qo) is 
that which minimizes the energy. It depends on 6. 

When the doping is very small Qq is close to tt, 7^2 =| cosQ |, and 



The minimum of the energy is at 



s 



t 



(5 + 0.079)2 J" ^^^^ 

The energy per site is plotted in Fig.l as a function of for different 
values of (5, t/J = 5 and 5 = 1/2. When 6 increases the minimum of the 
energy moves smoothly. A second minimum at Q = 27r appears when the 
doping is large enough. This minimum decreases when 6 increases and at 
some value of the doping the both minimums become equal. At this point Qo 
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has a step-wise behavior. The curves Qo{5) are plotted in Fig. 2 for different 
values oit/ J and S = 1/2. We treat the appearance of the second minimum 
as a appearance of instability and as a upper bound of applicability of our 
approach. For t/J= 10 this bound is (5 = 0.31, for t/ J = 5; 6 = 0.39, for 
t/J = S;6 = 0.47 and for t/J = 0, 5; J = 0.6. 

III. The spin-spin correlator 

Let us consider the spin-spin correlation function in the limit of large 
separation between the two sites. At zero temperature the bose system con- 
denses. Then, the large distance asymptotic of the bose correlators is deter- 
mined by the values of their Fourier transforms at the wave vectors at which 
the system condenses. From Eqs. (36-39) one obtains 

< d^aj >< a'j'tti > + < aittj >< a^Oj >~ (44) 

where 

\ 

m^ = -{m\ + ml) (45) 

The fermi term goes to zero when | — Vj \ goes to infinity and we shall 
drop it. As a result we obtain that the large distance asymptotic of spin-spin 
correlator at T = is 

< Si ■ Sj >= cos Qo ■ in - Tj) (46) 

This implies the presence of the long-range spiral phase. From Eqs. (26) 
one obtains that in leading order of 5 

mi = m2 = m = S -O.m -]-5 (47) 

In Fig. 3 we have plotted m as a function of 5 for different values oft/ J. 

At finite temperature, the bose energy Ei){k) (see Eq. (23)) has a finite 
gap. As a consequence, the bose correlators have a finite correlation length. 
For the spin-spin correlation function at small doping (when tt < Qo{S) < 
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Q2{5) and large distance one obtains 



< Si ■ Sj > 



1 cos Qo ■ {Vi - Vj) 



6e «i + 



{2SJoY 27r I n - Tj 




(48) 



(cos Qo) 



where 



ii' = 4V^i - 1 



(49) 
(50) 



-1 



4V(^2- I COsQo I)/ I COsQo 



From Eqs. (25) it follows that and ^2 grow exponentially when T goes 
to zero. At very small doping 7^1 ~ 7^2, ^1 — ^2, | cos Qo\—l and the spin-spin 
correlator looks like in the theory of the Heisenberg antiferromagnet^. When 
8 increases m2 decreases and the first term in Eq.(48) is dominated when T 
approaches zero. 

IV. Concluding remarks 

In conclusion, we have investigated the two-dimensional t — J model 
by modified spin-wave theory. A spin-spin correlator, we have obtained, at 
T = Eq. (46) shows a spiral long-range order. At T 7^ the correlation 
function decays exponentially as in the theory of the Heisenberg antiferro- 
magnet. The results, at small doping, are in qualitative agreement with the 
results obtained by Schwinger boson slave-fermion mean field theory^, except 
for the vector 3/2 known from the earlier publications. Substantial quanti- 
tative discrepancies appear when 5 increases. In spite of this the approach 
we have discussed above seams to be appropriate for the investigation of the 
normal state of the novel superconductors. 

As a by-product we obtain the spin-stiffness constant ps of the square 
lattice Heisenberg antiferromagnets. It can be defined by 



PS = 



Q=TT 
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where E{Q) is the energy per site of the system. From Eq. (42) aX 6 = one 
obtains 

PS = 2J{S + 0.079)^ 
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FIGURE CAPTIONS 



Fig. 1. The energy of the system per site as a function of Q for J = 5, 
5 = 1/2 and different values of doping. 

Fig. 2. The physical value Qo of the spirality angle as a function of 

doping, for t/ J = 10 (sohd line with circles), t/ J = 5 (dashed line), t/ J = S 
(dotted line), t/J = 0.5 (dot-dashed line) and S = 1/2. 

Fig. 3. The quantity m as a function of doping for t/J = 10 (solid line), 
t/J = 5 (dashed line), t/ J = 3 (dotted line), t/ J = 0.5 (dot-dashed line) and 
5 = 1/2. 
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